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A THREE-DIMENSIONAL ANALYSIS OF THE EFFECTS 


OF HEAT GENERATION IN PLATES 

By Robert R. McWithey 
Langley Research Center 

SUMMARY 

A method for determining three-dimensional steady- state temperature distributions 
is presented for plates in which heat is generated and conducted throughout the plate 
material. The solution for the temperature distribution is given in terms of a power 
series with respect to the plate thickness coordinate. The coefficients of the power 
series are dependent upon the temperature distribution over the plate surfaces. Specific 
attention is given to the determination of temperature distributions and plate deforma- 
tions resulting from constant heat generation and from the absorption of electromagnetic 
radiation. 


INTRODUCTION 

The phenomenon of heat generation in materials occurs as a result of the conver- 
sion of energy into heat in such processes as electrical-current conduction, chemical 
reactions, and absorption of electromagnetic radiation. Heat generation within a struc- 
tural material causes nonlinear temperature gradients to develop, and may result in 
undesirable thermal deformations or high thermal stresses. Some recent investigations 
have been carried out (refs. 1 to 6) to determine transient one-dimensional temperature 
distributions resulting from the absorption of various forms of electromagnetic radia- 
tion. In these analyses special attention was given to the determination of the thermal 
stresses and deformations in plate and shell configurations resulting from nuclear radia- 
tion (refs. 1 to 4). Examination of the equations governing plate and shell deformations, 
however, indicates that thermal deformation is also dependent upon the temperature 
gradients in the directions normal to the thickness coordinate. (See refs. 3 and 7 to 10.) 

It is therefore of interest to determine not only the one-dimensional (i.e., thickness direc- 
tion) but also the three-dimensional temperature distributions throughout the structural 
material. 

The work presented herein develops a method for determining three-dimensional 
steady-state temperature distributions caused by internal heat generation in plates. The 
corresponding expressions for the thermal-deformation parameters are presented for 



plates and shells. Specific examples are given for heat generation caused either by 
constant heat generation or by absorption of electromagnetic radiation. 

SYMBOLS 

a n (x,y),b n (x,y) coefficients of power series in z 

A(x,y,z) heat generated at any point within the plate per unit time per unit volume 
b 0 = T(x,y,0) 
bi = f (x,y,o) 

B n Bernoulli numbers 

Cn^rn constants 

E n Euler numbers 

f(x,y,z) = - Afepl 

h plate thickness 

I intensity of electromagnetic radiation 

I Q energy incident on plate surface z = h/2 

k thermal conductivity of plate material 

m,n integers in power series 

P o = |( T h/2 “ T -h/2) 

T temperature change from an initial state 

T h/2 = T(x,y,h/2> 

T _h/2 = T ( x »y,- h /2) 
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T ° " l( T h/2 + T -h/2) 


x,y,z 


rectangular coordinates 


(3 linear absorption coefficient 

fti = (2 2n - 2)B n 

a 


v i 2 


v 1 ” 

v n 

<P Q 

<p 1 


differential operator 
differential operator 


2 9 2 9 2 
8x 2 + 8y 2 + 8? 

82 q2 

8j? + 8y 2 


differential operator (v2) m / 2 ; V® indicates no differential operation 
intended 

differential operator (v2) ' • v® indicates no differential operation 
intended 

r h/2 

plate temperature parameter \ T dz 

J -h/2 

r h/2 

plate temperature parameter \ Tz dz 

J -h/2 


THEORY 


The following analysis presents a general series approach for determining three- 
dimensional steady-state temperature distribution in plates. The general series solu- 
tion is then used, in specific applications, for the determination of plate temperature dis- 
tributions and the corresponding thermal -deformation parameters. 

General Solution of the Heat- Conduction Equation 

Steady- state temperature distributions within solid bodies in which conduction is 
the mode of heat transfer and which have isotropic thermal properties are governed 
mathematically by Poisson's equation (ref. 11): 

v 1 2 t = - A fo y> ? )- (1) 

1 k 
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where 


V 1 




9 2 




T 


temperature at any point 


A(x,y,z) heat generated at any point within the plate per unit time per unit volume 
k thermal conductivity of plate material (assumed constant) 


Solutions of Poisson's equation may be found by means of the Green's function 
technique. (See ref. 12.) However, application of this method of solution to boundary- 
value problems is difficult since the solution involves a surface and volume integral over 
the body geometry. Another analytical method of approach to the solution of Poisson's 
equation may be formulated by expressing both the temperature and the right-hand side 
of equation (1) in terms of power series in z. This gives 


and 


OO 

T = £ b n (x,y)z n 

n=0, 1,2,3 


( 2 ) 


OO 

- A ( x k y,z ) = f(x,y,z) = ^ a n (x,y)z n 


(3) 


n=0 


where functions a n (x,y) are known because f(x,y,z) is assumed to be known. Substi- 
tuting equations (2) and (3) into equation (1) and equating the coefficients of like powers 
of z gives the recurrence relation 

a n = v2fe n + ( n + !)( n + 2)b n+2 (4) 


With the use of this recurrence relation, the coefficients, b n , of z n in equation (2) 
may be expressed in terms of bQ, bj, and a. n . Equation (2) then becomes 


- I <-*> 

n=0,2,4 


n/2 


A V n b 0 z n + 


,n+2 


11=0,2,4 


i <-» 

m=0,2,4 


m/2 (n - m)! v m 
(n+2): 


a n-m + | (-D (n ' 1 > /2 £v n ' 1 V n 

n=l,3,5 


* i * n * 2 S 

n=1 .3 R m=1 ! 


/ 1 \(m-l)/2 (n - m + 1): v m-l 
K ’ (n + 2).’ n-m+1 


(5) 
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where bg and bj are, respectively, the temperature distribution at the plane z = 0 
and the temperature gradient in the z-direction at the plane z = 0. If the functions bp 
and bj are known, the temperature distribution within the body is determined by equa- 
tion (5). The functions bq and bj are determined from the temperature boundary 
conditions. 


Determination of Coefficients bo and bj in the General Solution 

For the plate shown in figure 1, the temperature distribution within the plate is 
largely dependent on the temperature boundary conditions on the large surfaces (z = ±h/2) 

and is relatively unaffected by the temperature 
boundary conditions along the edges of the 
plate. Therefore, the edge boundary condi- 
tions will not be considered and the coeffi- 
cients bQ and bj will be functions of the 
boundary temperatures T ^/2 and T-h/2 
on the large surfaces of the plate and of the 
internal heating function a n (x,y). 



Figure L- Coordinate system associated with plate. 


In order to determine bo and bj in 
terms of Tj^* ^--h/2> an d an( x >y) it is 
convenient to define two boundary-temperature 
functions T 0 and P Q in terms of the speci- 
fied boundary temperatures 

( 6 ) 


T ° = 2 ( T h/2 + T -h/2) 

Po = ^ ("^h/ 2 - T -h/ 2 ) CO 

Expressing T Q and P Q in terms of bQ, bj_, and a n by the use of equation (5) gives 


T = T * + 
A o o + 


00 

1 (r 


n+2 


n=0,2,4 
n+2 


n 

V / i\ m /2 (n - m): v m 

L. } ’ (n + 2): a n 


m=0,2,4 


-m 


P 0 = Po + 


I (1 

n=l,3,5 


Y / 1 v(m-l)/2 (n - m + 1): „m-l 
L (n + 2): a n- 

m=l,3,5 


m+1 


( 8 ) 


( 9 ) 


where 


T 0 *» l (.l,»/2J ? v“b 0 (|) n 
n=0,2,4 N ' 


( 10 ) 
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and 


Po 


* 


£ (_l)(n-l)/2 

n=l,3,5 


_ 1 _ 

n! 



(ID 


It should be noted that equations (8) and (9) establish relationships between T 0 and bg 
independent of bj, and between P Q and bj independent of bg. Expressions for bg 
and bj may be obtained in terms of T 0 * and P 0 *, respectively, as outlined in 
appendix A. These expressions are given by 


b o= 2 

n=0,2,4 


( 12 ) 


n— 3, 5, 7 

where the constants E n are the Euler numbers, and the constants are re l a l e< l 

to the Bernoulli numbers (see eq. (A9)). 

Thus the functions bg and bj may be found in terms of T 0 , P 0 , and a n (x,y) 
by first determining T 0 * and P Q * from equations (8) and (9), and substituting the 
resulting expressions into equations (12) and (13). When the solutions for bg and bj 
are substituted into equation (5), the resulting temperature distribution will satisfy the 
governing Poisson equation (eq. (1)) and the specified temperature boundary conditions 
on the large surfaces of the plate. Appendix B presents general expressions for the 
temperature distribution in which Th/2> ^-h/2> anc * coefficients a n (x,y) in the 
series expansion for the heat- generation function are limited to either biharmonic func- 
tions or harmonic functions. 


Determination of the Plate Thermal-Deformation Parameters 

Equations governing the thermal deformations of plates may be found in numerous 
places in the literature. (See refs. 3 and 7 to 10.) The basic assumptions used in these 
analyses require the definition of two temperature parameters which appear in the equi- 
librium equations and boundary conditions governing plate deformation. These param- 
eters are given by 


c h/2 

(pn = \ T dz 

0 -h/2 

(14) 

C h/2 

w i = \ Tz dz 

1 -h/2 

(15) 
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and are proportional, respectively, to the average temperature distribution through the 
plate thickness and the first moment of the temperature distribution with respect to the 
midplane of the plate. For small-deflection plate theory the parameter <pq appears 
only in the equations for midplane stresses and in-plane deformations. Similarly, the 
parameter <p j appears only in the equations for bending stresses and out-of-plane 
deformations. (See ref. 10.) 

Expressions for the plate thermal-deformation parameters may be obtained in 
terms of bp and bj by substituting equation (5) into equations (14) and (15) and per- 
forming the indicated integrations. The resulting expressions for cpq and cpi are 
given by 


<Pq “ ^ 0 * + ^ 


I (M 

n=0,2,4 


,n+3 


V , ^m/2 (n - m)! „ 

L K ’ (n + 2): 
n + 3 m=0,2,4 


n-m 


(16) 


/h\ n+4 

<p, = cp* + 2 y ili— t (■i)(m-l)/2 (n-m + l): v m-l 

1 1 _ r'n c n + 4 A, r (n + 2): n-m+1 


(17) 


where 


and 


^0* = 2 J (-l) n/2 

n=0,2,4 


vnb o (hf +1 

(n + 1): \2 ) 


(18) 


<Pl 


* = 2 


(_ l )( n - !)/2 


v n-1 b 


n.— 1 ? 3 , 5 


(n + 2)(n!) 


a! ) \2/ 


,n+2 


(19) 


The functions cp^ and cp ^ may thus be determined from equations (16) and (17) after 
obtaining expressions for bo and bj which are functions of the temperature boundary 
conditions and the heating function. 

From the expression for the plate thermal-deformation parameter <pQ, it can be 
shown that (pQ is dependent only on the boundary temperature function T 0 and the 
coefficients a n for even values of n. Similarly, it can be shown that the plate thermal- 
deformation parameter cp ^ is dependent only on Po and the coefficients a n for odd 
values of n. Therefore, when the heating function A is represented by an even 
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function with respect to the midplane of the plate, the plate thermal-deformation param- 
eter is unaffected by heat generation. 

APPLICATIONS 


Constant Heat Generation 


General solution .- For the case of constant heat generation, the coefficients 
for the series defined by equation (3) are 


a =-^ 
a 0 k 

a n = 0 


(n > 0)( 


( 20 ) 


The temperature distribution given by equation (5) thus reduces to 


will „n oo I— n 1. n n 

. £ ( l) n/2 ° , £ ( ^(n-ll/2 v V Az 2 


n=0,2,4 


n=l,3,5 


n! 


2k 


( 21 ) 


The relationships between T 0 and bg and between P Q and b^ are obtained from 
equations (8) and (9) and are given by 


"\ 


rn rp * h A 

T °“ T ° "IT 


P 5= P 
*0 r o 


( 22 ) 


The thermal- deformation parameters are obtained from equations (16) and (17) and are 
given by 


= ,n * 


<P0 ~ <P o 
<P\ = <Pl 



(23) 


Equations (23) indicate that <Pq is the only plate temperature parameter affected by 
constant heat generation. Thus, for small deformations, since A is an even function, 
the constant internal heat generation affects only the in-plane plate deformation. In addi- 
tion, the effect of constant internal heat generation is to change the value of cpQ over 
the entire plate by a constant amount. 

2 

Harmonic solution .- If, in addition to a constant A, it is assumed that V T ^/2 

and V T_j ^2 are zero > t hen a simpler form of plate temperature distribution is 
obtained from equation (B16) of appendix B as 
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T =i(r- z2 ) + 


(24) 


2P q z 


+ T 0 


The corresponding plate thermal- deformation parameters are obtained from equa- 
tions (B17) and (B18) and are given by 


m u . h 3 A 

^0 o 12k 


and 


<P 1 = 


P h 
o 


(25) 

(26) 


Absorption of Electromagnetic Radiation 


General solution .- An assumption commonly made in calculations involving the 
absorption of electromagnetic radiation is that dl, the incremental amount of energy 
absorbed in an incremental thickness dz, is proportional to both the intensity of the 
electromagnetic radiation I and the thickness dz. (See ref. 13.) With the use of this 
assumption, the heat generated per unit time per unit volume is given by 


A(x,y,z) = 0I o (x,y)e 



(27) 


where I 0 is the energy incident and normal to the plane surface z = h/2 with no losses 
due to reflection and /3 is the linear absorption coefficient. 

The heating function A as described by equation (27) is zero for either j3 = 0, 

/3 — <», or I Q = 0. When equation (27) is substituted into equation (3), the summation 
of the power series in z with coefficients a n is defined as 



z n = 


n 


0I o (x,y) -0(|-z) 

k e 


(28) 


Expanding the right-hand side of equation (28) in a power series in z and equating like 
powers of z gives the expression for a n as 


an(x,y) = - 


Io<K.y>e-' Jh/ V +1 


knl 


(29) 


Substitution of equation (29) into equation (5) gives the temperature equation as 
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V n_ *bjZ n 


OU --yXX, uy T-TlA a, 

T = X <-l) n/2 ^ n + I (.1,(-1)/2L_5 
n=0,2,4 n=l,3,5 


■j8h/2 


00 

1 <-« 


n=0,2,4 


n/2 ^o 
/3n+l| 


,/3z 


n+1 

2 

m=0, 1,2,3 


m 


ml 


m! 


( 30 ) 


In order to satisfy the specified temperature boundary conditions, the functions bg 
and b^ in equation (30) must be determined in terms of T 0 * and P 0 * from equa- 
tions (12) and (13). First the functions To* and Po* are found by substituting equa- 
tion (29) into equations (8) and (9). After some algebraic manipulation, this gives T 0 * 
and P 0 * as 


t °*'T 0+ 4 2 


n/2 


n=0,2,4 


on+1 


e-f + X - 2e-^ 2 


V J -.M 

Lj m! V 2 / 

o a \ ' 


m=0,2,4 


(31) 


and 


OO 

P * = P +— Y (-l) n / 2 — 3 
P ° ° 2k L ( ’ gn+1 
n=0,2,4 p 


-e^ h + l-2e-^ 2 


n+1 


I ^l(f 

m=l,3,5 


m 


(32) 


The functions bg and bj can then be found by substituting equations (31) and (32) into 
equations (12) and (13). 

The plate thermal-deformation parameters may be found by appropriately inte- 
grating equation (30) for cpQ and cpi or by substituting equation (29) into equations (16) 
and (17). Either method gives 


<?0 = V 


2e 


-/3h/2 


OO 

2 <-» 


n/2 v %l 


0 n+2 


n=0,2,4 


sinh ^ - 


2. . (f) 


m=0,2,4 


m+1 


(m + 1): 


(33) 


cp 1 = C p 1 


2e 


■jSh/2 


OO 

2 <-» 


n/2 V % 


n=0,2,4 


0 n+3 


& cosh - sinh 
2 2 2 


n+1 

2 

m=l,3,5 


(m\ 

\ 2 / 


m+2 


(m + 2)ml 


(34) 
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Equations (33) and (34) are the general expressions for cpQ and cp^ in terms of bq 
and bj. (See eqs. (18) and (19).) Appendix B presents a solution which illustrates the 
general procedure used in obtaining T, q>Q, and <pj in terms of T 0 , P Q , and I 0 . 

Solution when ^’-h/2’ an< * *o are harmonic.- Two specific cases are pre- 

sented, one for temperature distributions in which T^/2 and ^-h/2 are prescribed, 
and one for temperature distributions in which Tjj/2 * s prescribed and the surface 
z * -h/2 is insulated. 


Temperature distributions in which T^/2 anc * ^-h/2 are prescribed: Appen- 
dix B presents the solutions for T, <Pq, and cp-^ for the case in which T^/2’ r *’-h/2> 
and I Q are biharmonic functions. (See eqs. (B19) to (B21).) If T^2> ^-h/2’ anc * *o 
are harmonic functions, the corresponding expressions for T, <pQ, and <p^ may be 
obtained from equations (B19) to (B21) by letting 


V I 2 T q = 0 
V 2 P 0 = 0 \ 

V 2 I o = 0 J 


(35) 


The resulting expressions are 


T = 




(36) 


I h 

<p 0 = T o h + — r~ 2 ^ h - 2 + e_ ^ h (^ h + 2)] 

u 2k(/3h)^ L J 


(37) 


cp 1 = 


Po h 


Ic)h — ^ |(j3h) 2 - 6/Sh + 12 - e"^ h [(/3h) 2 + 6/3h + lJl 
12k(/3h) 3 l J 


(38) 


Equations (37) and (38) indicate that the plate temperature parameters are linearly 
dependent on the temperature distributions over the large surfaces of the plate and on 
the magnitude of the incident radiation. 

The value of the linear absorption coefficient j3 also affects the magnitude of the 
plate thermal-deformation parameters inasmuch as it governs the value of the coeffi- 
cient of the incident- radiation term. Figures 2(a) and 2(b) show the effect of /3h on 
coefficients of the incident- radiation terms from equations (37) and (38). It may be seen 
from figure 2(a), where (cp^ - T 0 h)kyi 0 h2 is plotted against /3h, that radiation absorp- 
tion has the largest effect on <pQ, and hence on the in-plane small-deflection plate 
deformations, for a value of /3h of approximately 2.6. It should be noted that cp^ is 


11 



3h 


(a) Effect of ph on ^g. 

Figure 2 Effect of ph on plate temperature parameters for specified boundary temperatures 
and incident energy distribution. 


shows the quantity cp 


Po h2 \ 


equ l al to T 0 h when /3h is equal to zero. Also, cp^ approaches T 0 h as /3h 
approaches infinity. In these limiting cases no radiation is absorbed within the plate, 
and the temperature gradient is constant through the thickness of the plate. Figure 2(b) 

jly^I 0 h 2 plotted against /3h, and indicates that the quan- 
tity which determines the out-of-plane small-deflection plate deformations, 
increases from a value of P 0 h^/6 for a value of /3h equal to zero until 0h is 7.14. 
Thereafter, decreases with increasing /3h and approaches P 0 hy/ 6 as 0h 
approaches infinity. As before, the extreme values of 0h are cases in which no radia- 
tion is absorbed within the plate and the temperature gradient is constant through the 
thickness of the plate. 


Temperature distributions in which T^/2 is prescribed and the surface z = -h/2 
is insulated: For this case the function T_ h /2 is eliminated from the T 0 and P Q 

terms of equation (36) by using the condition 


8T| 

9z z= -h/2 


= 0 


(39) 
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The resulting expression for the temperature distribution is 


T = T 


h/2 




,-ffo 


Z - 


, 4 -§) 


/3 




( 40 ) 


Because no heat flow is allowed across the surface z = -h/2, the temperature gradient 
through the thickness of the plate must allow heat conduction toward the surface z = h/2. 
Thus the maximum temperature in the thickness direction always occurs at the surface 
z = -h/2. A typical temperature profile, as obtained from equation (40), is shown in 
figure 3 in terms of the dimensionless temperature parameter ^T - Tj^j^Ao 11 anc * the 

dimensionless thickness parameter 2z/h. The magnitude of the temperature difference 
between the front and back surfaces of the plate is determined by the parameter /3h. 


Figure 4 shows the quantity ^T ~ ^h/2) V* 0 * 1 Pitted as a function of /3h 
from equation (40). It is seen that a maximum temperature difference between the front 
and back surfaces will occur at a value of /3h equal to 1.8. Thus, appropriate selec- 
tion of plate material and plate thickness will allow temperatures at the back surface 
(z = -h/2) to be maximized with respect to the front surface temperature. 



-1.0 - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1.0 

2_z_ 

h 


Figure 3.- Temperature profile through plate thickness when jp= 0 at z = j) and ph = 1.0. 
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The plate temperature parameters <?q and cp^ are obtained by appropriately 
integrating equation (40). The resulting equations are 

O’o ' h V2 + { e ^ h t 2 - <® 2 ] * ^ - 2 } < 41 > 

= J o h3 [e'^lm 3 - 60h - 12] + 1 - ft (42) 

1 k(/3h)3^12 2 j 

Figure 5(a) shows the variation of (cp^ - hT^g)^/!^^ with /3h as obtained from equa- 
tion (41), and indicates that the largest effect of radiation absorption on the value of cp q 
occurs for /3h - 2. Figure 5(b) shows the variation of <pjk/i 0 h3 with /3h as obtained 
from equation (42), and indicates that the maximum effect of radiation absorption on the 
value of occurs when the value of 0h is 1.61. 


15 



(a) Effect of ph on <p Q . 

Figure 5.~ Effect of ph on the plate temperature parameters when the 
distribution are specified on the surface z = h/2, and wf 
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temperature distribution and incident energy 
len dT/dz = 0 at z = -h/2. 



0 2 4 6 8 10 12 14 16 18 20 

Ph 

(b) Effect of ph on <p\. 

Figure 5.- Concluded. 

CONCLUDING REMARKS 

General expressions for the steady-state temperature distribution within plates 
have been obtained for conditions involving heat generation within the plate. The solu- 
tions presented satisfy the temperature boundary conditions on the large surfaces of the 
plate but do not usually satisfy the temperature boundary conditions at the edges of the 
plate. 

The two temperature parameters that govern plate deformation are denoted by cp^ 
and cp^ and are respectively proportional to the average temperature through the plate 
thickness and the first moment of the temperature distribution in the thickness direction 
with respect to the midplane of the plate. The equations indicate that when the internal 
heating function is represented by an even function with respect to the midplane of the 
plate, the parameter cp j is unaffected by heat generation. The solution for constant 
heat generation is presented as an example of this case. If the internal heating function 
is represented by an odd function, </>q is unaffected by heat generation. 
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When the generation of heat is induced by the absorption of electromagnetic radia- 
tion, the resulting nonlinear temperature distribution through the plate thickness affects 
both the plate deformation parameters. The variation of the deformation parameters 
with respect to the dimensionless product /3h (/3 is the linear absorption coefficient 
and h is the plate thickness) is of particular interest inasmuch as the magnitudes of /3 
and h are dependent on the plate material and plate thickness, respectively. For the 
case in which the temperature boundary conditions and the incident radiation are har- 
monic functions, it is observed that 0h has the greatest effect on the magnitudes of <?q 
and (p i for values of /3h between zero and approximately 7. Therefore, if it is 
desirable to minimize the effect of internal absorption of electromagnetic radiation on 
thermal deformation, it is necessary to select a material and thickness for which the 
magnitude of the product /3h is either very nearly zero or much greater than 7. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 28, 1967, 

124 - 08 - 01 - 13 - 23 . 
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APPENDIX A 


DETERMINATION OF COEFFICIENTS bg AND b 1 

The temperature distribution within the plate as defined by equation (5) is dependent 
on the unknown midplane temperature functions bg and bj. Equations (10) and (11) 
relate bg and bj to the known boundary temperature functions and heat- generating 
function but are not directly solvable for bg and b^ unless V^bg = V^bj = 0. Expres- 
sions for bg and bj are found in terms of the known functions T 0 * and P Q *, how- 
ever, by appropriately reversing the series equations for T 0 * and P Q *. For example, 
a set of equations involving T 0 * and its derivatives may be obtained in the following 
form: 



where C n are arbitrary constants. If these constants satisfy the set of equations 
given by 


C 2 - 1 = 0 
C 2 (4)(3) 


Ct . 

„ C 4 (6)(5) . Qb(6)(5) , n 

c 6-^r- + ^r“- 1 = 0 

„ c 6 (8)(7) C 4 (8) (7) (6) (5) C 2 (8)(7) , A 

C 8 22 + 4! 2 — +1 -° 


(A2) 




19 


APPENDIX A 


then the summation of equations (Al) is simply 



(A3) 


Equations (A2) are of the form presented in reference 14, series (1145), and it is 
immediately apparent that equations (A2) are satisfied if 

^n = ®n (A4) 

where the constants E n are Euler numbers and are defined in reference 14, pages 238 
and 239. Thus, equation (A3) becomes 


bQ 


- 1 % 
n=0,2,4 


, ,n 

h^ v n T 


(A5) 


In determining bj in terms of P 0 , a set of equations may be written in the form 


, * bih v \/h \ 3 + v S/h \ 5 

0 ~ 2 3! V2/ 5! \2y 



+ 


V 8 b, 


9: 





Dl/h' 

2 :\ 2 , 



* 


v*b 


l/h\3 


v 4 b 


1 2: U; “ D i 2:3: I2 


+ D 


v \hf D 

1 2:5: 12) 1 


Aw 9 . 

2:7: 12J 


£2/h 

4: 12 




+ d 2 


v 8 b. 


4: 5: 



> (A6) 
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where the constants D m are arbitrary constants. If these constants satisfy the set of 
equations given by 
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Dl 

— - 0 

2: 

3: 

1 

CM | 

Q 1 

£L + A- 

4 : 

2:3: 5: 

D3 

D 2 Dj 

6: 

4:3: 2 : 5 : 

d 4 

d 3 D 2 

8: 

6:3: 4 : 5 : 


i- = 0 


D 1 1 
■=-=- + -=-= 0 


(A7) 


then the summation of equations (Al) is simply 

= P * + Y D ( n -l)/2 ftf- 1 (n-l), ♦ 

2 0 Jx, <” - W 

By using series (1135) of reference 14 along with the identities 

/3 n = (2 2n - 2) B n for all n 
B n (l) = 0 for all n 

( 1 ' 


(A8) 


B n ±1=0 


for odd n 

J 


(A9) 


which are listed respectively in series (1130), (1136), and (1142) of reference 14, it may 
be shown that, for odd values of n!3 


0 = - 4 + 


h 


$2 


n! 2I(n-2): 41 (n - 4): 


+ • • . 


/ i\(n- 1)/2 q 
- (-1) P(n-l)/2 


(n- 1)! 


(A10) 


Equations obtained from (A10) have the same form as equations (A7) . Equations (A7) 
are therefore satisfied if 

Vl)/2 S Vl)/2 (All) 

where the constants are re l a ^ e< ^ to the Bernoulli numbers as indicated by 

equation (A9). Thus equation (A8) becomes 


b, = 


2P f 


y ^(n-i)/2 /h\ n -y-i r 

k n ^ W 


(A12) 


n — 3, 5,7 
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SIMPLIFICATION OF TEMPERATURE EQUATIONS 


Temperature Distributions in Which T h y 2 > T -h/2’ and ^ 

Are Biharmonic Functions 

When completely arbitrary temperature boundary conditions and heating function 
are introduced into the general series solution for the temperature distribution (eq. (5)) 
and into the series for b 0 and bj (eqs. (12) and (13)), the mechanics of solution 
become unwieldy. However, if the boundary temperatures T h /2> T -h/2> and the 
heating-function coefficients a n (x,y) are restricted so that 

v mr 


L h/2 

-h/2 


V mr p 


V m a n = 0 



(Bl) 


where m is some reasonably small, even integer, then the solution for the temperature 


distribution is more easily obtainable. For example, if it is assumed that T^/2) T -h/2> 

and a n are biharmonic functions then, 


v4T h/2 - 
V 4 T. h/2 = 0 / 

(B2) 


< 

tP 

II 

O 

(B3) 

Applying equations (B2) and (B3) to equations (6) to (9), (12), and (13) gives: 


v 4 t 0 = 0 

(B4) 


V 4 P 0 = o 

(B5) 


v 4 t 0 * = 0 

(B6) 


O 

11 

* 

O 

t> 

(B7) 


v 4 b 0 = o 

(B8) 


V 4 b : = 0 

(B9) 

Substituting equations (B3), (B8), and (B9) into the general equation for the temperature 
distribution (eq. (5)) gives 

T = b 0 - V ^°z2 + b]Z - V 6 b l z 3 + ^z2 + *l z 3 


oo 

+ Y z n+2 

y2a (n-2) 

(BIO) 

+ Z 

n=2,3,4 

(n + 2)(n +1) (n + 2)(n + l)(n)(n - 1) 
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where bQ and b ^ are obtained from equations (12) and (13) as 


(Bll) 


2P n h 2 * 

b i = -fr + ii vp o 


(B12) 


Expressions for T 0 and P Q are obtained from equations (8) and (9) in terms of T 0 , 
P Q , and a n and are given by 

,2 ^ /, \n+2 


T„ = T n - 


aph _ 

8 n=2^4,6^"^ 


oo 

1 .81 




V 2 a 


(n-2) 


(n + 2)(n +1) (n + 2)(n + l)(n)(n - 1) 


(B13) 


3-1 h /Vi\] 

°* = p ° " ~lo~ " X J2 ) 


n=3,5,7 


\n+2 


n 


V 2 a 


i Sl2L 


(n + 2)(n + 1) (n + 2)(n + l)(n)(n - 1) 


(B14) 


The temperature distribution may now be written in terms of T 0 , P 0 , and a n by 
eliminating T 0 * and P 0 * from the expressions for bg and bj (eqs. (Bll) and 
(B12)) and substituting the resulting expressions for bg and bj into equation (BIO). 
This procedure gives the temperature distribution as 

V / h\ n+2 a n V /h\ n+4 v2& n 

T " T ° ' n= ^ 2 4 W (n + 2)(n + 1) + n= Z 2 4 \2/ (n + 4)(n + 3)(n + 2)(n + 1) 


( I ) 2 


v 2 t 


o_i y /h\ n * 2 , 2 r7 2 r 

2 n= o 2 4 ' (n + 2)(n+l) 3h 


_z_ V /h\ n+2 v > 

3h /-! \2 ) (n + 2)(n + 1) 

n=l,3,5 ' ' 


OO 

♦ i - 

n=0,l,2 


n+2 an 


00 

1 


.n+4 


v2a n 


n=0,l,2 


+ — 

(n + 4)(n + 3)(n + 2)(n +1) h 


(n + 2)(n + 1) 


Po- I (If 

[_ n=l,3,5 x ' 


an 


V /h\ n+4 v2a n 

n=l 3 5 X ' ( n + 4 )( n + 3 ^ n + 2 H n + !) 


(n + 2)(n + 1) 

(B15) 
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Expressions for <p^ and <p^ may be obtained by appropriately integrating equa- 
tion (B15) with respect to z (see eqs. (14) and (15)). 


Solutions in Which T^g* ^’-h/2 , an< * % Are Harmonic Functions 

If the further restriction is made that T'-h/2» and a n are harmonic func- 

tions, then the v2 terms appearing in equations (BIO) to (B15) vanish and the tempera- 
ture distribution becomes 


2P n z 

T = T 0 + — — - + 
u h 


y *nh 2 - (in , . y a * 1 - or 

A, . (n + l)(n + 2) n=l 3 5 ( n + 1 ^ n + 2 ) 


n=0,2,4 


Applying these results to the plate thermal deformation parameters <pQ and cp^ 
(eqs. (14) and (15)) gives 

v (i 

^0 c T ° h Z 4 (nVl)(n,3) 

n=0,2,4 


= T o 


sr 4 *. 


V.2 y 1,2; ^ 

^■*•6 3 Z/ (n + 2)(n + 4) 

n=l,3,5 


Biharmonic Solution for the Case of Absorption of Electromagnetic Radiation 

The temperature distributions caused by the absorption of electromagnetic radia- 
tion are obtained by substituting the expression for a n , given by equation (29), into 
equation (B15). This gives 



where it should be noted that V^i 0 must be zero because of the requirement that V^a n 
be zero (see eq. (29)). Appropriately integrating equation (B19) with respect to z gives 
(pQ and < p^ as 
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